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ABSTRACT 

We analyze the effect of dissipation on the orbital evolution of supermassive black holes (SMBHs) using high- 
resolution self-consistent gasdynamical simulations of binary equal- and unequal-mass mergers of disk galaxies. 
The galaxy models are consistent with the ACDM paradigm of structure formation and the simulations include 
the effects of radiative cooling and star formation. We find that equal-mass mergers always lead to the formation 
of a close SMBH pair at the center of the remnant with separations limited solely by the adopted force resolution 
of ~ 100 pc. Instead, the final SMBH separation in unequal-mass mergers depends sensitively on how the central 
structure of the merging galaxies is modified by dissipation. In the absence of dissipation, the satellite galaxy can 
be entirely disrupted before the merger is completed leaving its SMBH wandering at a distance too far from the 
center of the remnant for the formation of a close pair. In contrast, we show that gas cooling facilitates the pairing 
process by increasing the resilience of the companion galaxy to tidal disruption. Moreover, we demonstrate that 
merging disk galaxies constructed to obey the Mbh-c relation, move relative to it depending on whether they 
undergo a dissipational or collisionless merger, regardless of the mass ratio of the merging systems. Collisionless 
simulations reveal that remnants tend to move away from the mean relation highlighting the role of gas-poor merg- 
ers as a possible source of scatter. In dissipational mergers, the interplay between strong gas inflows associated 
with the formation of massive nuclear disks and the consumption of gas by star formation provides the necessary 
fuel to the SMBHs and allows the merger remnants to satisfy the relation. 



Subject headings: black hole physics 
numerical 

1. INTRODUCTION 



cosmology: theory — galaxies: mergers — hydrodynamics — methods: 



Irrefutable dynamical evidence indicates that supermassive 
black holes (SMBHs) with masses ranging from 10 6 to above 
10 9 M Q reside at the centers of most galaxies hosting spheroids 
(e.g., Kormendy & Richstone 1995). The available data show 
the existence of a remarkably tight correlation between the mass 
of the SMBHs, Mbh, and the stellar velocity dispersion of the 
host galaxy spheroidal component, a (Ferrarese & Merritt 2000; 
Gebhardt et al. 2000), suggesting a fundamental mechanism 
that connects SMBH assembly and galaxy formation (e.g., Silk 
& Rees 1998; Burkert & Silk 2001). 

According to the currently favored cold dark matter cosmo- 
logical models, structures in the universe are the end result of 
a complex hierarchy of mergers and accretion of smaller sub- 
units. Thus, the hierarchical buildup of SMBHs by massive 
seed black holes present at the center of protogalaxies and the 
formation of SMBH binaries appear as a natural consequence 
in any "bottom-up" cosmogony. Recently, the discovery of a bi- 
nary active galactic nucleus in the interacting system NGC 6240 
by Komossa et al. (2003) has lent support to this picture. Al- 
though substantial gas accretion is also required for the SMBHs 
to reach their present-day mass density and satisfy the Mbh-c 
relation ( Volonteri et al. 2003), the link between these two main 
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modes of SMBH growth is still missing. 

The formation and subsequent orbital evolution of a SMBH 
pair depends on how efficiently the host galactic cores lose 
angular momentum by dynamical friction during the merging 
process. Once the two cores have merged, the SMBHs may 
decay further and eventually form a binary owing to the drag 
exerted by the background mass distribution. These two evolu- 
tionary phases have been investigated numerically by a number 
of authors (Governato et al. 1994; Makino & Ebisuzaki 1996; 
Milosavljevic & Merritt 2001; Makino & Funato 2004). How- 
ever, the galaxy models used in these studies were idealized, 
spherical stellar systems that could at most faithfully represent 
the central spheroidal components of real galaxies. Hence, both 
the larger scale dynamical evolution of the merging systems and 
the cosmological framework were missing. Moreover, these 
studies did not explore the role of a dissipative component in 
driving the evolution of a SMBH pair. Notable exceptions are 
the studies by Escala et al. (2004a,b), who showed that gas 
causes continuing loss of angular momentum and rapidly re- 
duces the relative SMBH separation to distances where gravi- 
tational radiation and coalescence is efficient. However, if and 
how their initial conditions are related to the larger scale dy- 
namics involved in galaxy merging is still unclear. 

In this Letter we report on the effects of gaseous dissipation 
on the fate of SMBHs using high-resolution binary disk galaxy 
merger simulations. Governato et al. (1994) have already high- 
lighted the difficulty of forming a close SMBH pair when tidal 
disruption of one of the host systems intervenes. As we illus- 
trate below, gaseous dissipation is vital for the survival of the 
host galaxies by deepening their potential wells and provides 
the necessary fuel for the SMBH growth. Finally, we analyze 
the stellar kinematics and distribution of gas in the central re- 
gions of the merger remnants to investigate for the first time 
how merging galaxies move with respect to their initial loca- 
tion on the Mbh-c plane. 
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FIG. 1 . — Spherically averaged density distribution for different components in the initial models (left) and some of the remnants (right) plotted from the force 
resolution outward. Left: Different thick curves show the total (solid), dark matter density after the adiabatic contraction due to baryons (dashed) and total density 
of baryons (dot-dashed) in the basic galaxy model. The dotted curve shows the dark matter density before the baryonic contraction. The thin solid curve shows the 
total density of the satellite galaxy. After the adiabatic contraction the central density slope is close to isothermal with p(r) oc f~ 2 for both galaxy models. Right: 
Total density profiles for some of the merger remnants in the coplanar encounter geometry. The numbers in brackets indicate the values for the gas fraction, / g , used 
in the particular simulation. The center of each remnant was defined as the location of the minimum of the gravitational potential. In dissipative simulations, the 
remnants exhibit an increase in the central density resulting from strong gas inflows and star formation associated with the merger. 



2. NUMERICAL SIMULATIONS 

We perform four types of high-resolution simulations of bi- 
nary disk galaxy mergers with mass ratios of 1 : 1 and 4: 1 includ- 
ing different physical processes. We simulate purely collision- 
less mergers and mergers in which we follow the gas dynam- 
ics "adiabatically," i.e., without radiative cooling. A third set 
of simulations includes radiative cooling, while in the fourth 
set we also allow the cold gas to form stars. The simulations 
were performed with GASOLINE, a multi-stepping, parallel 
TreeSPH JV-body code (Wadsley et al. 2004). We include radia- 
tive and Compton cooling for a primordial mixture of hydrogen 
and helium. The star formation algorithm is based on that by 
Katz (1992), where gas particles in dense, cold Jeans unstable 
regions and in convergent flows spawn star particles at a rate 
proportional to the local dynamical time (see also Governato 
et al. 2004), and reproduces the Schmidt law. The star forma- 
tion efficiency was set to 0. 1 which yields a star formation rate 
of 1 -2 M Q yr _1 for models in isolation that have a disk gas mass 
and surface density comparable to those of the Milky Way. 

The multicomponent galaxy models are constructed using 
the technique described by Hernquist (1993) and their structural 
parameters are consistent with the ACDM paradigm for struc- 
ture formation (Mo et al. 1998). To this end, each galaxy con- 
sists of a spherical and isotropic Navarro et al. (1996) dark mat- 
ter (DM) halo (Kazantzidis et al. 2004a), an exponential disk, 
and a spherical, Hernquist (1990) non-rotating bulge. For the 
basic galaxy model we adopted parameters from the Milky Way 
model Al of Klypin et al. (2002). Specifically, the DM halo 
had a virial mass of M v ; r = 10 12 M Q , a concentration parame- 
ter of c = 12, and a dimensionless spin parameter of A = 0.031. 
The mass, thickness and resulting scale length of the disk were 
M d = 0.04M vir , zq = 0.1R A , and R d = 3.5 kpc, respectively. The 
bulge mass and scale radius were Mb = 0.008M v ; r and a = Q.2R A , 
respectively. The DM halo was adiabatically contracted to re- 
spond to the growth of the disk and bulge (Blumenthal et al. 
1986) resulting into galaxy models with a central total density 
slope close to isothermal (Figure 1). The companion galaxy 



is either a replica of the same model in equal-mass mergers 
or a system containing one-fourth of the mass in each compo- 
nent in which lengths and velocities are renormalized accord- 
ing to the cosmological scaling for virialized systems (Mo et al. 
1998). The disk scale length of the satellite galaxy is equal to 
Rd = 2.2 kpc. Given the uncertainties in the adopted M/L, our 
galaxies are consistent with the stellar mass Tully-Fisher and 
size-mass relations. To each of the galaxy models we added a 
particle representing a SMBH at the center of the bulge compo- 
nent that was initially at rest. For the larger galaxy model we 
used a SMBH mass equal to M BH = 3 x 10 6 M Q . 

We considered two values for the gas fraction, / g , namely 
10% and 50% of the total disk mass. We shut off radiative 
cooling at temperatures below 2 x 10 4 K that is about a fac- 
tor of 2 higher than the temperature at which atomic radiative 
cooling would drop sharply due to the adopted cooling func- 
tion. With this choice we take into account non-thermal pres- 
sures and approximate the simulated gas with the warm ISM 
of a real galaxy (Barnes 2002). Moreover, this choice enables 
us to simulate very gas-rich disks that would otherwise become 
strongly gravitationally unstable and undergo widespread frag- 
mentation. Two main encounter geometries were adopted: pro- 
grade coplanar mergers and mergers in which one of the disks 
was inclined with respect to the orbital plane. We also simu- 
lated a single retrograde equal-mass merger in which both disk 
spin vectors were antiparallel to the orbital angular momentum 
vector. The galaxies approached each other on parabolic orbits 
with pericentric distances that were 20% of the more massive 
galaxy's virial radius, typical of cosmological mergers (Khoch- 
far & Burkert 2003). The initial separation of the halo centers 
was twice their virial radii and their initial relative velocity was 
determined from the corresponding Keplerian orbit of two point 
masses. In collisionless mergers, we used N = 1.2 x 10 6 parti- 
cles in total, with each galaxy consisting of 10 5 disk particles, 
10 5 bulge particles, and 10 6 DM particles. The gas component 
in gasdynamical simulations was represented by 10 5 particles. 
We adopted a gravitational softening of e = 0.1 kpc for both 
the DM and baryonic particles of the larger galaxy, and for its 
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FIG. 2. — Final separation of SMBHs (filled circles) in a subset of merger 
simulations. The large scale (left) and small-scale (right) structure of the rem- 
nants projected onto the orbital plane is also shown. All frames correspond 
to remnants that were allowed to relax for several dynamical times after the 
merger was complete. The top and bottom rows of panels present results for 
the coplanar 1:1 and 4:1 mergers, respectively, with gas cooling. The mid- 
dle rows of panels corresponds to the coplanar collisionless 4: 1 merger. The 
frames on the left show the logarithmic baryonic surface density maps and are 
320 X 230 kpc. The limiting surface density is S = IMq/ pc 2 . Blue and red 
maps are used for the stellar and gaseous component, respectively, and adap- 
tive smoothing is used to preserve details in high-density regions. The top and 
bottom frames on the right show the central gaseous disks, while the middle 
frame on the right shows the stellar distribution of the merger remnant. 



SMBH e = 0.03 kpc, which is small enough to follow its orbital 
evolution. For the satellite galaxy, all softening lengths were 
scaled according to e oc m 1 / 3 . In all simulations the merger rem- 
nants were allowed to settle into equilibrium ~ 15 dynamical 
times at the disk half-mass radius after the merger was com- 
pleted. 

3. RESULTS 

3.1. SMBH Pairing in Binary Disk Galaxy Mergers 

The galaxies merge in three to five orbits (between 5.5 to 7 
Gyr) depending on the mass of the companion, a timescale that 
is set by dynamical friction and tidal stripping (e.g., Taffoni 
et al. 2003). These timescales are longer than those reported in 
earlier studies (e.g., Barnes & Hernquist 1996) due to the larger 
and more realistic pericentric distances adopted here. 

In equal-mass mergers, the cuspy potentials of both galaxies 
are deep enough to allow the survival of their inner regions un- 
til orbital decay by dynamical friction is complete. This result 
holds for both collisionless and gasdynamical mergers. In the 
latter, during the first two orbits a strong spiral pattern appears 
in both the stellar and the gaseous component and mild non- 
axisymmetric torques redistribute mass and angular momentum 
driving approximately 20% of the gas towards the center. How- 
ever, the central bulges stabilize the galaxies against bar forma- 
tion. Shortly before the galaxies merge, a second much stronger 
gas inflow occurs caused by strong tidal torquing as the galax- 
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ies perform a close fly-by. However, only in the simulations 
with radiative cooling does this inflow result in a significant 
central concentration of cold (r~2x 10 4 K) gas. In this case, 
more than 80% of the gas that was originally in the disks is col- 
lected within the central 500 pc. This results in a considerable 
steepening of the central density slope of the combined distri- 
bution of DM and baryons (Figure 1). In adiabatic simulations, 
only 20-30% of the initial gas mass collects within the central 
kpc. This gas is heated to temperatures T > 10 s K by shocks 
and compressions during the merger and subsequently expands 
since it cannot radiate away its energy. 

Large central gas inflows were already noticed in lower reso- 
lution merger simulations (e.g. Barnes & Hernquist 1996; Barnes 
2002). In simulations with star formation these large inflows re- 
sult in a central starburst and more than 90% of the central gas 
distribution is converted into stars in less than 10 8 yr. The peak 
star formation rates range from 30 to > 100 M Q yr _I , compara- 
ble to those of luminous infrared galaxies (LIRGs) and ultralu- 
minous infrared galaxies (ULIRGs), suggesting that our mod- 
eling of star formation and the amount of cold gas present are 
realistic. The two SMBHs end up orbiting at the center of the 
remnant on eccentric orbits and form a close pair at a separation 
comparable to the adopted force resolution of ~ 100 pc (Fig- 
ure 2). Significantly higher mass and force resolution would be 
needed to follow the SMBHs at relative separations where they 
would form a binary system. 

In unequal-mass mergers, the outcome depends sensitively 
on how the internal structure of the merging galaxies is mod- 
ified by dissipation. In collisionless simulations, the tidal dis- 
ruption of the satellite galaxy at about 6 kpc from the center 
leaves its SMBH wandering at a distance that prohibits the for- 
mation of a close pair. This "naked" SMBH would contribute 
to a population of wandering SMBHs in galactic halos (Volon- 
teri et al. 2003) At such a distance the timescale for orbital 
decay predicted by dynamical friction would exceed a Hubble 
time. In simulations with cooling, the situation is completely 
reversed as the gas inflow becomes particularly strong in the 
companion galaxy. We measured an inflow of about 3 M Q yr _1 
within the central ~ kpc compared with < 0.5 M Q yr _1 for the 
larger galaxy. This striking difference in the gas inflows is at- 
tributed to a tidally induced strong stellar bar observed in the 
satellite galaxy approximately 1 .5 Gyr before the final passage. 
More than 50% of the gas is funneled into the center owing to 
gas shocking in the bar potential. At this stage the companion 
galaxy has about 20% more gas within its central ~ kpc rel- 
ative to the larger galaxy. Dissipation enables the core of the 
infalling galaxy to survive complete tidal disruption by deep- 
ening its potential well resulting in a SMBH pair separated by 
~ 100 pc as in equal-mass mergers (Figure 2). 

The large inflows observed in the cooling and star forma- 
tion simulations always produce a rotationally supported nu- 
clear disk with a size in the range 1-2 kpc and a tempera- 
ture set by our temperature floor (T = 2 x 10 4 K). In adiabatic 
simulations, this gas remains in a hot phase and forms a pres- 
sure supported cloud of similar size. The nuclear disks result 
from the coalescence of two central disk-like structures formed 
at the centers of the galaxies primarily by the strong inflow 
just before they merge. These disks are tilted by several to a 
few tens of degrees relative to the orbital plane of the galax- 
ies. They have peak rotational velocities and gas masses in the 
range of 250 -300 km s" 1 and ~ 10 s - 10 9 M Q , respectively, and 
they are well-resolved (A^IO 4 ) in our high mass resolution sim- 
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ulations. Nuclear disks of molecular gas of a few hundred pc 
to a ~ kpc scale have been identified spectroscopically for sev- 
eral AGNs and ULIRGs (e.g. Downes & Solomon 1998; Davies 
et al. 2004). The observed sizes, masses, and rotational veloci- 
ties are comparable to those measured in our simulations. The 
nuclear disks reported here provide the reservoir of gas that fu- 
els the central SMBHs. 

3.2. Mbh-c Relation in Binary Disk Galaxy Mergers 

For each remnant, we calculate line-of-sight aperture stellar 
velocity dispersions a within one effective radius R e (Gebhardt 
et al. 2000) and we use a range of viewing orientations, namely 
inclinations i = 0, 60°, and 90°. For i = 60° and i = 90°, we 
rotate the system about an axis perpendicular to the inclination 
axis through an angle <f> set to 0, 45°, 90°, and 135°. Thus, 
for each simulation we obtain 9 different orientations. For each 
of these viewing angles, we measured the mass density profiles 
in circular apertures out to two radii (8 and 25 kpc) under two 
assumptions: (1) the system is a disk galaxy and decompose 
the mass density profile into a Sersic bulge and an exponential 
disk; (2) the system is a pure spheroid and fit it with a pure 
Sersic profile. This procedure ensures that the measurements 
of R e , and therefore of a, that we obtain are very robust. 

In cases when the two SMBHs form a close pair, the Mbh 
assigned to each remnant is calculated by summing the masses 
of the two SMBHs and adding the gas mass contained within 
the resolution of our gasdynamical simulations. In doing this 
we implicitly assume that any close pair will actually form a 
SMBH binary that will subsequently coalesce. Indeed, this 
is supported by the smaller scale simulations of Escala et al. 
(2004a,b), who showed that SMBH pairs starting at a separa- 
tion of few hundred pc in a gaseous background will typically 
merge within ~ 10 7 yr. When the SMBHs do not form a close 
pair, we only consider the mass of the SMBH located closer to 
the center of the remnant. Figure 3 shows the location of all 
merger remnants on the Mbh-c plane together with real galaxy 
data and their best-fit correlation taken from Tremaine et al. 
(2002). We note that the initial galaxy models start very close 
to the mean correlation by construction. The error bars show the 
spread about the mean value of a and are smaller in the cooling 
and star formation simulations likely due to the fact that these 
remnants are more spherical (Kazantzidis et al. 2004b). Merg- 
ers lead to an increase of the stellar velocity dispersion that is 
remarkably different depending on the physics included in the 
simulations. The largest increase occurs in simulations with ra- 
diative cooling which exhibit the deepest potentials. When star 
formation is included, we measure velocity dispersions closer 
to that of the collisionless simulations, since the gas does not 
build a significant central mass concentration. 

Figure 3 illustrates that collisionless mergers move remnants 
away from the mean Mbh-c relation. In adiabatic simulations, 
the remnants have in principle enough gaseous fuel to remain 
close to the relation, but the gas is too hot to accrete onto the 
SMBHs. In mergers with radiative cooling, the total amount 
of nuclear gas is more than one order of magnitude larger than 
needed to keep the galaxies close to the relation. This should 
also be regarded as an upper limit on the gas mass that may ac- 
crete onto the SMBHs since feedback processes ensuing once 
the AGN becomes active will regulate the inflow and allow the 
accretion of only a fraction of this gas, thus limiting the SMBH 
growth (Silk & Rees 1998). If AGN feedback is strong enough 
to supress cooling a hot gas phase similar to that found in our 
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FIG. 3. — Mbh-c relation in binary disk galaxy mergers. The open triangles 
show data from the galaxy sample compiled by Tremaine et al. (2002) and the 
solid line corresponds to their best-fit correlation. Results for simulated rem- 
nants are shown in color and the initial galaxy models are denoted by stars. 
Open and filled symbols correspond to collisionless and gasdynamical simula- 
tions, respectively. Circles and squares show equal- and unequal-mass merg- 
ers, respectively. The filled triangle corresponds to an unequal-mass merger 
in which the larger disk was inclined by 45° with respect to the orbital plane. 
Symbols outlined with crosses correspond to mergers with 50% gas fraction in 
the disks. The error bars show the spread about the mean value of a in each 
remnant. 



adiabatic simulations is expected to form. Thus, the results 
from cooling and adiabatic simulations likely provide upper and 
lower bounds to the true physical behavior of the resulting rem- 
nants. Note that the point corresponding to the 4: 1 merger is the 
highest with respect to the best-fit correlation in Figure 3. This 
is the product of the particularly strong gas inflow occurring in 
the satellite galaxy just before the merger is completed. Overall 
a higher fraction of the available gas is driven towards the cen- 
ter relative to the equal-mass galaxy mergers. When star for- 
mation is included the amount of cold gas in the nuclear disks 
is reduced by more than 90% and the merger remnants move 
nearly parallel to the relation. This suggests that star formation 
during mergers is a key ingredient for maintaining the tightness 
of the Mbh-c relation. 

4. DISCUSSION AND CONCLUSIONS 

Gaseous dissipation influences considerably the outcome of 
binary mergers of disk galaxies containing SMBHs. Most im- 
portantly, it controls the SMBH pairing process in unequal- 
mass mergers by modifying the central structure of the com- 
panion galaxy, enabling it to survive complete tidal disruption. 
Dissipationless, unequal-mass mergers are also expected to pro- 
duce close SMBH pairs when the satellite galaxy is dominated 
by a dense stellar central component, as in nucleated dwarf el- 
liptical (e.g. M32) and SO galaxies. Our results suggest that 
semi-analytic models of hierarchical SMBH growth that ne- 
glect the effect of dissipation on the orbital evolution of SMBHs 
underestimate their pairing efficiency and subsequent coales- 
cence. This consequently leads into overestimating the number 
of wandering SMBHs in MW-sized galaxies (Volonteri et al. 
2003). The higher pairing efficiency of - 1O 6 M SMBHs re- 
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ported here has important implications for the probability of ob- 
serving coalescence events with space-based gravitational waves 
experiments such as LIS A(Sesana et al. 2004). 

Gaseous dissipation also determines how merging galaxies 
constructed to satisfy the M Bli -<r relation move relative to their 
initial location on the M Bli -cr plane. The collisionless simula- 
tions reveal that when galaxies become gas poor which is likely 
at low z, their merger remnants will tend to move away from the 
mean relation. This suggests that gas-poor mergers act as a pos- 
sible source of scatter in the mean Mbh-c relation. The natural 
scaling between the amount of cold nuclear gas and the increase 
in the central stellar velocity dispersion explains why merger 
remnants in simulations with radiative cooling move above the 
mean relation. Star formation serves to move merger remnants 
nearly parallel to the relation. The fact that unequal-mass merg- 
ers appear effective at building a reservoir of gas for SMBH ac- 
cretion is noteworthy, since in hierarchical structure formation 
models they are significantly more frequent than equal-mass 
ones. 

The nuclear disks at the center of the remnants show signifi- 
cant spiral patterns (Figure 2) which transport angular momen- 
tum outwards and sustain radial inflows towards the center. We 
measure typical inflows in the range 10~ 2 - 10~ 3 M Q yr _I within 
1 kpc. The spiral instabilities are sustained by the self-gravity of 
the disks. In star formation simulations which contain the light- 
est disks, the minimum Toomre Q parameter (Toomre 1964) 
is ~ 2 within the same distance. Since gravity is softened at 
scales corresponding to a significant fraction of the disk size 
(200 pc), short wavelength modes are stabilized and the disks 
behave effectively as if they had a higher Q parameter (Mayer 
et al. 2004). Therefore, spiral modes should be more vigor- 
ous with increasing resolution and promote even stronger in- 
flows. Such inflows could feed the SMBHs and bridge the gap 
between the large scale flows and the viscous accretion taking 
place once the gas has reached the AGN accretion disk at sub- 
parsec scales (Heller & Shlosman 1994; Fukuda et al. 2000). 
Future simulations will be used to explore the feasibility of such 
mechanism. Interestingly, observed nuclear disks also have Q 
values between 1 and 2 and some of them exhibit significant 
non-axisymmetric structure (Downes & Solomon 1998). 
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